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In the preceding paper, we assembled the theoretical components necessary for a unified framework of singlet fission, 
a type of multiexciton generation producing two triplet excitons from one singlet exciton. In this paper, we apply our 
methodology to molecular dimers of pentacene, a widely studied material that exhibits singlet fission. We address 
a longstanding theoretical issue, namely whether singlet fission proceeds via two sequential electron transfer steps 
mediated by a charge-transfer state or via a direct two-electron transfer process. We find evidence for a superexchange 
mediated mechanism, whereby the fission process proceeds through virtual charge-transfer states which may be very 
high in energy. In particular, this mechanism predicts efficient singlet fission on the sub-picosecond timescale, in 
reasonable agreement with experiment. We investigate the role played by molecular vibrations in mediating relaxation 
and decoherence, finding that different physically reasonable forms for the bath relaxation function give similar results. 
We also examine the competing direct coupling mechanism and find it to yield fission rates slower in comparison 
with the superexchange mechanism for the dimer. We discuss implications for crystalline pentacene, including the 
limitations of the dimer model. 



I. INTRODUCTION 

Singlet exciton fission, a process whereby one high-energy 
singlet exciton is converted into two lower energy triplet exci- 
tons, is an excited state phenomenon with a potential impact 
on the efficiency of inexpensive organic solar cells. 1 In light 
of this potential utility in the design of photovoltaic systems, 
there has been a recent explosion of experimental studies on a 
variety of molecular materials 2-10 and fabricated devices 
Currently, a detailed microscopic understanding of this pro- 
cess is lacking. In the preceding paper, 14 we outlined a fully 
microscopic theoretical framework for the practical simula- 
tion of singlet fission dynamics. In particular, within the con- 
text of singlet fission chromophore systems, we connected ex- 
cited state quantum chemistry with established reduced den- 
sity matrix methods from quantum relaxation theory. 

Here, we continue this endeavor with a realistic theoretical 
treatment of the singlet fission dynamics of molecular dimers. 
Understanding singlet fission in small molecular complexes 
has intrinsic interest and such complexes may have utility as 
a sensitizer in a Gratzel-type solar cell^i^ They may also 
be representative model systems to understand the singlet fis- 
sion process in bulk crystals. We focus on pentacene because 
it is perhaps the most thoroughly studied material that has 
been robustly shown to exhibit singlet fission. In particular, 
experimental evidence points to fast and efficient singlet fis- 
sion in bulk pentacene, 5 *^ - — but the process appears to be 
much slower in the molecular (tetracene) complexes realized 
so far. 20,21 Here we will explore dimer systems both as a real- 
istic description of small molecular complexes and as a model 
potentially applicable in bulk materials. Full application of 
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the present formalism to pentacene clusters and crystals will 
be subject of a future report. 

An ongoing debate, the crux of which is laid out in Ref. 1 
and in our previous paper^ concerns the quantum mechani- 
cal mechanism by which singlet fission takes place amongst 
the manifold of excited states. Specifically, this debate rests 
on the role, or lack thereof, of charge-transfer (CT) states in 
the singlet fission process. The relevant quantum mechani- 
cal states as well as the two competing mechanisms are de- 
picted schematically in Fig. [TJ The so-called "mediated" 
mechanism, posits that a single electron transfer from the in- 
tramolecular excited singlet state, ISiSo) or l^o^iX produces 
a charge-transfer state, \CA) or \AC), after which a second 
one-electron transfer event yields the multi-exciton, triplet- 
triplet state, \T\T\). Alternatively, a "direct" mechanism im- 
plies a simultaneous two-electron process which circumvents 
CT states altogether. 




FIG. 1. Schematic representation of the five electronic states relevant for 
singlet fission in a dimer. The actual states employed in the calculations are 
spin-adapted linear combinations yielding overall spin-singlets, unlike those 
shown here. 
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Previous theoretical work on singlet fission has largely con- 
sidered only the mediated mechanism. In particular, Greyson 
et al. investigated the case of purely coherent energy transfer 
for mediated fission 22 and Teichen and Eaves derived solvent- 
dependent rate expressions for the separate one-electron trans- 
fer events implicit in the mediated mechanism. 23 The latter 
authors concluded that a necessary condition for efficient me- 
diated fission is that the CT state energy lie in between that of 
the intramolecular singlet and the multiexciton triplet-triplet 
state, i.e. E(SiS ),E(S S{) > E(CA),E(AC) > E(T X T X ). 

This viewpoint is consistent with a recent quantum chem- 
istry calculation on pentacene clusters reported by Zimmer- 
man et al.*2i They concluded that because CT states were cal- 
culated to be significantly higher in energy than intramolecu- 
lar singlets (by about 300 meV or more), singlet fission in pen- 
tacene cannot take place via the mediated mechanism. Rather, 
the authors supported the direct mechanism, estimating a di- 
rect coupling matrix element of about 5 meV. However, this 
number is almost two orders of magnitude smaller than that 
required to explain the experimentally observed timescale of 
fission in pentacene, ranging from 80 - 200 fs i 5 A 17 ~ — 

The present manuscript is a first step towards the resolution 
of this apparent paradox. Because singlet fission is inherently 
a dynamical process, one must exercise caution in the inter- 
pretation of static electronic structure calculations and their 
implications for fission. As such, we argue that a microscopic, 
dynamical treatment of the relevant electronic states coupled 
to a finite temperature bath is crucial for a theoretically sound 
description of singlet fission processes. Furthermore, the ac- 
curacy of the methodology and its associated approximations 
must be established for these complex problems. The imple- 
mentation should be carefully benchmarked and thoughtfully 
parametrized for the relevant physical problem, in this case, 
singlet fission. We have carried out the first step of benchmark 
calculations in our previous paper and here we take the second 
step, parametrizing a system-bath Hamiltonian for fission in 
molecular dimers and using an accurate quantum relaxation 
master equation to calculate the fission dynamics. Through 
this program, we are able to make firm statements regarding 
the feasibility of competing mechanisms as well as predict and 
rationalize experimental fission rates. 

The layout of the paper is as follows. We begin in Sec. HU 
with a review of the methodology presented in our previous 
paper. 14 In Sec. [Till we present our results for pentacene, which 
explore the effects of energy levels, electronic couplings, and 
phonon properties. Although we use pentacene as an example 
molecule, our exploration of important singlet fission param- 
eters is sufficiently broad so as to elucidate generic aspects 
of singlet fission. We summarize our work and conclude in 
Sec. En 



II. METHODOLOGY 

In this section, we briefly describe the adopted theoretical 
methodology as laid out in our previous paper, to which the 
reader is referred for more details.— In essence, we employ a 
system-bath Hamiltonian describing the coupling of the elec- 



tron and phonon degrees of freedom, 25,26 H tot = H e i + H e i- P h + 
ft ph , with 



H el = 2 imai + 2 ii>vyoi, 



(1) 



Hei-ph = ^ 1001 ^ Ck,m + ^ 1001 ^ c kJj q k , (2) 



and 



H ph = ^ 



1 



Pk , 1 2-2 

y + 2^* 



(3) 



The parameters of this Hamiltonian are determined via a vari- 
ety of ab initio and semi-empirical methods, and the dynamics 
generated under the action of this Hamiltonian are calculated 
by a perturbative quantum master equation. The electronic 
structure and quantum dynamics methodologies are described 
in the following sections. 



A. Geometry and electronic structure 

To evaluate the role of molecular geometry, we consider 
individual pentacene dimers extracted from the ab plane 
of the experimental crystal structure. 27 There are approxi- 
mately three symmetry unique nearest-neighbor dimer pairs, 
as shown in Fig.O corresponding to the [a b] translation vec- 
tors [1 0], [1/2 1/2], and [-1/2 1/2]. By investigating these 
dimers in particular, we are probing the extent to which effi- 
cient singlet fission observed in crystals but not dimers may 
be due simply to molecular orientation as opposed to ex- 
plicit many-body environmental effects. Alternatively, for 
a different choice of electronic structure parameters (which 
may be altered by electronic polarization for example), the 
dimers may be taken as model systems for the crystal, as- 
suming purely local energy transfer events as described in 




FIG. 2. Molecular geometry of the pentacene crystal. Three pentacene 
molecules are emphasized, displaying the three symmetry-unique nearest- 
neighbor dimer pairs discussed in the text. Also shown are isosurface plots 
of the HF HOMO (a) and LUMO (b) of the isolated molecules, including the 
phase convention adopted in this work. 
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the introduction. Though such a picture permeates the lit- 
erature, one should not invest too much in this model-based 
view, instead preferring a direct treatment of the full system 
as done in our forthcoming work on pentacene crystals. Exci- 
ton derealization in neat acene crystals may span up to tens 
of molecules, 8 invalidating this simple dimer picture. Before 
concluding, we will discuss analogous calculations on a cova- 
lently linked pentacene dimer inspired by previously studied 
tetracene dimer a 2Q i 21 with the aim of rationalizing the low ob- 
served fission yield. 

When applied to a molecular dimer, the electronic model 
Hamiltonian described in our previous paper yields five dia- 
batic states (excluding the ground state), each one a spin sin- 
glet. The first two are localized, intramolecular Frenkel exci- 
tations, |SiSo) an d l^o^i)- Additional single excitations gen- 
erate the third and fourth states of charge-transfer character, 
\CA) and \AC). The fifth and final state is a doubly excited 
triplet- triplet state \T\T\), presumed to be the spin singlet pre- 
cursor to fully separated triplets. These five states are depicted 
schematically in Fig.ffl though we emphasize that the single- 
configuration states shown there are not states of well-defined 
spin multiplicity, i.e. they are not eigenfunctions of the S 2 
operator. The spin-adapted variant of these dimer states em- 
ployed here (as well as their non- spin- adapted counterparts) 
have been used as a starting point in a variety of other theoret- 
ical works . l i 22 i 23 

Because calculated excited state energies exhibit errors on 
the order of 0.5 eV, especially for acenes, 28 we take the diag- 
onal matrix elements of the electronic Hamiltonian, E[, to be 
adjustable parameters. This procedure avoids any bias inher- 
ited from electronic structure methodology and more impor- 
tantly provides for qualitative, physical insight into the effect 
that electronic energies have on singlet fission dynamics. The 
electronic couplings, Vy, on the other hand, are taken from 
ab initio calculations using the Hartree-Fock (HF) molecu- 
lar orbitals (MOs) of isolated pentacene molecules. As laid 
out in App. A of our previous paper and Eqs. (7)-(13) of 
Ref. 01 one-electron couplings are given by off-diagonal el- 
ements of the Fock operator for the combined, two-molecule 
system. Specifically, we have 

(CA\H e i\SiS ) = t LL , (AC\H e i\S 1 S ) = -t HH , 

{CA\H e i\T\T\) = V3/2 t LH , {AC\H el \T x T x ) = ^Jl t HL , 

(CA\H d \SoSi) = -t H H, (AC\H d \S Si) = t LL . 

In the above, t LL and t HH denote the one-electron coupling 
of the LUMO and HOMO, respectively, whereas t LH is the 
electronic coupling between the LUMO of the first molecule 
(A) and the HOMO of the second (B) and likewise for tHL, i.e. 



t L L = (L A \F\L B \ 
t LH = (L A \F\H B ), 



t HH = (H A \F\H B ), 
t HL = (H A \F\L B ). 



It should be noted that the mixed couplings, tuj and tHL, are 
entirely responsible for CT-mediated singlet fission. Thus, 
whereas previous theoretical studies of singlet fission have in- 
voked Longuet-Higgins-type approximations to estimate the 
mixed couplings, 22 the direct evaluation in terms of the Fock 
operator employed here should be preferred. 



Additionally, the direct evaluation yields the sign of all 
electronic couplings, in contrast to approximate methods 
(such as the energy- splitting in a dimer method 29 ) that yield 
only the magnitude of the coupling (however the sign can of- 
ten be inferred by inspection of the orbitals). In multistate 
systems and band theory calculations, the sign of the coupling 
can be very important and should be retained whenever possi- 
ble. To achieve a consistent sign, one must adopt a phase con- 
vention for the molecular orbitals of the system. In agreement 
with previous studies^ we employ an approximate screw- 
axis to fix the phase of the MOs, which are shown as used 
in Fig. [2 

In addition to the one-electron coupling matrix elements, a 
variety of two-electron integrals are apparent in the coupling 
expressions given in our previous paper and Ref. 1. These 
two-electron integrals, when included, were calculated by rep- 
resenting the HF MOs on a real-space grid and subsequently 
utilizing fast Fourier transform techniques. 31 Integrals were 
checked for convergence with respect to the grid size. 

All quantum chemistry calculations were performed with 
either the GAMESS (US) quantum chemistry package 32 or the 
Firefly quantum chemistry package, 33 which is partially based 
on the GAMESS (US) source code. Calculations employed 
the 6-31G(d) basis set. 



B. Quantum dynamics 

Reduced density matrix (RDM) quantum dynamics calcula- 
tions were performed within the Redfield framework ^ 25 i 26 i 34 ~ — 
with the secular and Markov approximations, whose use was 
justified theoretically and numerically in our previous paper 14 
to which we refer the reader for technical details. Briefly, 
the electronic RDM, which follows from a trace of the to- 
tal density matrix, W(t), over the phonon degrees of freedom, 
Pij(t) = (i\Tr p hW(t)\j), describes the time evolution of the pop- 
ulations, pa = Pi, and coherences, pij, of the electronic states. 
In the adiabatic basis that diagonalizes the electronic Hamilto- 
nian, H e i\a) = hcL) a \a), the RDM obeys the equation of motion 



dp a p(t) . Vl 

— - -ICOapPapyt) + 2_j K af3ySpyd(t). 

y,6 



(4) 



The Redfield tensor, 5 35 R a p 7 d, which effects finite- 
temperature relaxation and dephasing processes, can be ex- 
pressed in terms of thermal bath correlation functions. Those 
correlation functions are entirely determined by the spectral 
densities of the phonon degrees of freedom, 



4M = ^ y — sic* - co k ) 

2 r 



and 



Jij(<j)) = - > 6(a) - o) k ), 

2 r <* 



(5) 



(6) 



which may be calculated through a combination of molecu- 
lar dynamics and quantum chemistry calculations. Thus, the 
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[a b] 


tHH 


t L L 


t H L 


tin 


[1 0] 


85 (34,51) 


-60 (-43) 


-74 


74 


[1/2 1/2] 


-145 (47,-74) 


116 (-82) 


109 


-124 


[-1/2 1/2] 


228 (-85,131) 


-111 (84) 


108 


-134 



TABLE I. Electronic coupling parameters (in meV) of pentacene for the three 
dimer types described in the text. Values in parentheses are those calculated 
by Yamagata et al^- and Troisi and Orlandi, 37 the latter only where available 

(tHH). 



form of the spectral density adopted in a RDM calculation en- 
capsulates the complicated detailed motion of the vibrational 
phonon degrees of freedom. For simplicity, most of our re- 
sults will employ the common Ohmic spectral density with a 
Lorentzian high-frequency cutoff, 

J°{a>) = 2^,.^— L- (7) 
J CO 1 + £2f 

I J 

with the the strength of the system-bath interaction quanti- 
fied by the reorganization energies Ay = n~ l J dojJ^(oj)/oj, 
and the frequency of the interaction quantified by the cutoff 
Qij. To investigate the crucial role played by phonons in our 
dynamical calculations, we will also consider two alternative 
forms of the spectral density, to be discussed in Sec. MI Cl and 
shown in Fig.[9j This latter study highlights a particularly im- 
portant advantage of the Redfield formalism, because it can 
easily be applied to any form of the spectral density, unlike 
many other methods. 

III. RESULTS FOR PENTACENE 

We now proceed to apply the above methodology to pen- 
tacene. The one-electron coupling parameters calculated as 
described are given in Tab.H and are in reasonable agreement 
with values obtained by semiempirical calculations 30 and den- 
sity functional theory (DFT). 37 All electronic coupling val- 
ues are clearly on the 100 meV order of magnitude, though 
we point out that HF appears to systematically yield larger 
couplings as compared to other methods. A uniform rescal- 
ing may be performed as in Sec. IIIIBi but we will not do so 
here. Henceforth, we will only present results for the [1/2 
1/2] dimer, but results for the other two are qualitatively sim- 
ilar, with quantitatively different dynamical timescales. As 
discussed above, the diagonal energies will be varied in our 
simulations, the results of which are presented in Sec. IIII A[ 
Furthermore, we initially neglect all two-electron integrals, 
thereby investigating only the mediated fission mechanism. 
Combining the computed one-electron couplings with the 
aforementioned expressions thus yields the electronic Hamil- 
tonian in units of meV, 
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116 


-152 


145 





116 


E(SiS ) 








145 


-152 
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145 
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145 


133 


116 
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The phonon bath will be characterized by only diagonal 
system-bath coupling, initially of the Ohmic form, Eq. (|7]), 
with Jfrjico) = 0, i.e. we exclude the off-diagonal Peierls 
coupling, but will revisit this topic in Sec. IIII Cl We have per- 
formed the calculations presented here in the presence of off- 
diagonal coupling to a low-frequency bath and find the results 
to be largely unchanged; such low-frequency modes are sim- 
ply not efficient at mediating large-scale energy transfer, i.e. 
fitly «c \Ei - Ej\. For simplicity, we assume identical, uncor- 
rected baths for each electronic state, i.e. £lu = Q and An = A. 
In reality, the bath parameters will be slightly different for 
each state and perhaps correlated because they share certain 
molecules. While we surmise that these effects will only be 
of quantitative significance, they are a potentially interesting 
topic for further research. The cutoff frequency is fi£l =180 
meV = 1450 cm -1 , as is typical for acenes and other conju- 
gated organic molecules. The reorganization energy is A = 50 
meV = 400 cm -1 . Though these parameters imply a Huang- 
Rhys factor, S = A/H£l ^0.3, which is about a factor of two 
lower than the experimental one (see our previous paper for a 
theoretical comparison to the experimental vibronic spectrum 
of Ref. 38), they insure a quantitatively accurate treatment of 
dynamics within the weak-coupling, Redfield framework, as 
demonstrated numerically in our previous paper. In Sec. IIII Cl 
we will consider larger, more realistic values of the reorgani- 
zation energy, for which the Redfield theory predictions are 
still expected to be qualitatively accurate. 

The initial condition of the reduced density matrix is p(0) = 

IS1S0XS1S0I, ^' e - on ty one m °l ecu l e is initially excited. Al- 
ternative initial conditions, such as a coherent or incoherent 
mixture of |SiSo> an d |SoSi>, were found to yield qualita- 
tively similar results. In most simulations, we find that this S \ 
superposition is achieved, regardless of the initial condition, 
on a 10 fs timescale, which is much faster than the overall 
fission process. 



A. Scanning energies 

Because we only consider a homodimer and exclude the 
ground state, we can simplify our notation, collectively re- 
ferring to the S i So and SoSi states as "Si", CA and AC as 
"CT\ and T\T\ as "7T". In this case, there are only two 
independent energetic parameters, which we take to be the 
energy offset of the CT states with respect to the TT state, 
E(CT) - E(TT), and the analogous offset of the S i states to 
the TT state, E(S i) - E(TT). We naturally expect that fission 
will take place as long as both of these parameters are positive, 
such that TT is the lowest energy state (an assumption well- 
founded for pentacene, based on experiment and calculations). 
This expectation is clearly validated in Fig. [51 which shows the 
singlet fission yield after the four periods of time, t = 0.1, 0.2, 
0.5, and 1 ps. The fission yield at time t is calculated simply as 
the population of the diabatic TT state, PttQ), times 200%, 
the latter factor indicating conversion of one exciton into two. 
Note that this metric is a combination of the rate of fission as 
well as the thermodynamic equilibrium of the system. One 
can easily imagine situations where the rate of fission may be 
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(a) Yield after* = 0.1 ps 



(b) Yield after t = 0.2 ps 




0.0 0.5 



(c) Yield after t = 0.5 ps 



0.0 0.5 1.0 1.5 
(d) Yield after t = 1 ps 




0.0 0.5 1.0 1.5 
E(Si) - E(TT) [eV] 



0.0 0.5 1.0 1.5 
E{S X )-E{TT) [eV] 



FIG. 3. Singlet fission yield, Prrit) x 200%, after the four periods of time 
indicated for the [1/2 1/2] pentacene dimer. The dashed line qualitatively 
separates the superexchange (SX) regime, Eqt > E$ 1 , from the sequential 
(SEQ) regime, Es l > Eqt . Estimated energy levels for the pentacene dimer 
are denoted by the white circle. 



fast but thermodynamic equilibrium does not overwhelmingly 
favor the TT state. Different physical situations and techno- 
logical applications will dictate whether it is more desirable 
to extract fewer carriers due to fission at short times (wherein 
one would want to optimize the rate only) or wait longer to 
extract more carriers (optimize the equilibrium). Of course 
the combination, i.e. rapid fission with high thermodynamic 
efficiency, is most ideal and may also be possible in some sit- 
uations. 

Returning to Fig.[3j we divide the energetic phase space into 
two regions, E(CT) >E(S\) and E(CT) < E(S\), demarcated 
by a dashed white line. Only the latter yields the rather obvi- 
ous energetic pathway for mediated fission, i.e. population 
flows from S i to CT, and then from CT to TT. We call this 
the "sequential" mediated mechanism. The sequential mech- 
anism can be clearly seen in Fig. [3] to yield very efficient sin- 
glet fission, even at short time. Perhaps the ideal energetic 
configuration is achieved for E(S i) - E(TT) = 400 meV and 
E(CT) - E{TT) = 200 meV, which yields about 150% singlet 
fission after only 200 fs. 

We now consider the opposite energetic regime, E(CT) > 
E(S\). Although this regime naively suggests a barrier to sin- 
glet fission (recall that the direct coupling term has been set 
to zero), we see a remarkably high fission yield, even at short 
times, as long as < E(Si) - E(TT) < 500 meV. With this 
criterion satisfied, efficient singlet fission occurs even for CT 
energies up to 1 eV above TT. We refer to this somewhat sur- 
prising result as the "superexchange" mediated mechanism, a 
phenomenon familiar from electron transfer in magnetic and 
photosynthetic systems, 39-41 and introduced in our preceding 
paper 14 in the context of singlet fission. Clearly near the 



boundary E(CT) = £(£1), the distinction between "sequen- 
tial" and "superexchange" is not so sharp. However we will 
continue to adopt these names, so as to imply that the dy- 
namics are mostly characteristic of either one or the other, i.e. 
these limiting forms provide a useful language for the discus- 
sion of competing effects in CT-mediated singlet fission. 

One may naturally question the relevance of the above anal- 
ysis to pentacene, asking what are the relevant energetic pa- 
rameters for a pentacene dimer? In all panels of Fig. [3l 
we have placed a circle that encompasses the estimated en- 
ergy levels for pentacene (discussed below), clearly placing 
it in the superexchange-dominated regime. Recent calcula- 
tions on pentacene dimers by Greyson et al* 2 ^ using a com- 
bination of time-dependent DFT and constrained DFT, found 
E(Si)- E(TT) = 240 meV and E(CT) - E(TT) = 354 meV. 
As pointed out, accurate electronic structure calculations of 
excited states, including those with multiple excitations, can 
be difficult and so we also consider estimates based on ex- 
perimental measurements. To a first approximation (which 
was also adopted by Greyson et al.), the energy of the multi- 
exciton state is simply twice the energy of the lowest triplet 
state, E(TT) « 2E(T X ) = 2 x 0.86 eV = 1.72 eV, where 
0.86 eV is the experimental T\ excitation energy. 42 The first 
singlet excitation energy of a pentacene monomer is approx- 
imately E(S\) = 2.1-2.3 eV 43,44 giving an energy offset of 
E(S\) - E(TT) « 400 meV. Diabatic charge transfer ener- 
gies are difficult to determine experimentally, but estimates 
from (adiabatic) spectroscopic measurements on crystals sug- 
gest values upwards of 2.3-2.5 eV, 30,45 thereby predicting 
E{CT) - E(TT) « 600 meV or more. 

Most importantly, there is little debate that charge transfer 
energies are always higher in energy than those of the first ex- 
cited singlet, such that pentacene lies unambiguously in the 
superexchange regime of CT-mediated singlet fission. While 
we of course cannot definitively conclude that fission in pen- 
tacene occurs exclusively via CT-mediated superexchange, 
(given the many approximations in our work and because we 
have not yet addressed the possibility of direct fission) we can 
say with certainty that high-lying CT intermediate states do 
not preclude efficient CT-mediated fission. 

To visually summarize the results of this section, we first 
show in Fig. |4] an energy level diagram depicting the mixing 
of diabatic states to form adiabatic states in the energetic ar- 
rangement E(S 1) - E(TT) = 250 meV and E(CT) - E(TT) = 
500 meV which is approximately correct for pentacene and 
characteristic of the superexchange regime. The calculated 
population dynamics for this system are shown in Fig. [5] in 
both the diabatic (a) and adiabatic (c) basis; the S\ and CT 
populations are given by Ps l s ( t )^ p s s l (0 and PcaW+PacQ) 
respectively. In the diabatic basis one observes a very-short 
time mixing of S 1 and CT, after which CT remains approx- 
imately constant while S 1 decays into TT with a single rate 
constant. This behavior is exactly that of conventional su- 
perexchange, although the CT population is slightly larger 
than typical due to the strong electronic coupling. The same 
behavior can be observed perhaps more directly in the adi- 
abatic basis, where the S\ <-> CT superposition is popula- 
tion near-instantaneously, which then decays to an adiabatic 
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diabatic adiabatic 

A S^CT 

2xcr= 



(a) Superexchange, diabatic (b) Sequential, diabatic 



250 meV 
2*S X = 
250 meV 
TT- 



162 meV 

—S^CT 

158 meV 
— TT 



FIG. 4. Energy level diagram depicting the diabatic electronic states (i.e. be- 
fore mixing) and the adiabatic electronic states (i.e. after mixing), for a typ- 
ical "superexchange" energy configuration indicative of a pentacene dimer. 
For adiabatic states which are a significant mixture of two different types of 
diabatic states, the notation i <-> j is employed. 



state of essentially TT character. These dynamics should be 
contrasted with those of a sequential fission mechanism, with 
E(Si)- E(TT) = 500 meV and E(CT) - E(TT) = 250 meV, 
shown in Fig. [2b) and (d). In both bases it is clear that a two- 
step kinetics prevails whereby an initially excited state first 
decays into an intermediate of CT character, which then it- 
self decays into the final TT state. Though both mechanisms 
yield highly efficient singlet fission on the 1 ps timescale, their 
underlying mechanistic details are clearly quite distinct. 



B. Superexchange and the strength of the electronic 
coupling 

To understand how superexchange arises, consider the first- 
order effect that coupling to CT states has upon the initially 
excited S \ states, 



\sf } ) + 



Vs u CT 



-icr (0) > 



E(Si)-E(CT) 
from which the effective coupling from S i to TT follows as 

Vs u ctVct,tt 



(9) 



(S\ l) \H el \TT^) 



E(Si)-E(CT) 

Vs u ctVct,tt 



[E(S - E(TT)] - 

-Vs u ctVct,tt 
[E(CT) - E(TT)] ' 



[E(CT) - E(TT)] 



(10) 



In the last line above we have assumed that E(CT)-E(TT) » 
E(S\) - E(TT). For the electronic parameters of pentacene 
considered above, this effective coupling is approximately 
50 meV. Performing second-order semiclassical (Marcus-like) 
perturbation theory in this effective electronic coupling yields 



k(Si TT) 



2n 



^/4nh 2 Ak B T [E(CT) - E(TT)f 
[E(S i) - E{TT) + A] 2 



(11) 
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FIG. 5. Population dynamics contrasting superexchange and sequential CT- 
mediated singlet fission, shown in both the diabatic and adiabatic representa- 
tions. Diabatic energy levels for panels (a) and (c) are £(Si) - E(TT) = 250 
meV, E(CT) - E(TT) = 500 meV; and for panels (b) and (d) are reversed, i.e. 
E(S i) - E(TT) = 500 meV, E(CT) - E(TT) = 250 meV. 



from which one can directly read off a sharp, Gaussian de- 
pendence on the S \ energy gap, with width « ^4Ak B T ^75 
meV, and a very weak, power-law dependence on the CT en- 
ergy gap, both of which are consistent with the data previously 
shown in Fig. [3] Furthermore, it is apparent that the superex- 
change mechanism is a process which is overall fourth-order 
in the electronic coupling matrix elements, Vy. 

In light of this analysis, one should expect that an alterna- 
tive master equation that is perturbative to second-order in the 
electronic coupling (akin to Forster theory or the noninteract- 
ing blip approximation (NIB A), as in the work of Teichen and 
Eaves in Ref . [23|) would only predict sequential, and not su- 
perexchange, CT-mediated fission mechanisms. In Fig.[6l we 
show the fission dynamics predicted by the NIBA-like master 
equation, 



dPjjt) 
dt 



dsK„(t, s)Pj(s), 



(12) 



where 



Knit, S) = 



n 2 



-Re (exp (-Mft/h) exp (iHffs/h)) 



ph 



and Hjf = (i\H tot \i), clearly demonstrating that such master 
equations incorrectly predict no fission if E(CT) > E(S i). 

To check the validity of the derived fourth-order scaling, 
we introduce a dimensionless parameter, r/, which character- 
izes the electronic coupling strength. Specifically, we replace 
Vij — > TjVij, and consider the limit r] — > 0, for which a su- 
perexchange mechanism predicts k oc if. In Fig. 13a), we 
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(a) Yield after t = 0. 1 ps (b) Yield after t = 0.2 ps 




(c) Yield after t = 0.5 ps (d) Yield after t = 1 ps 




E(St) - E(TT) [eV] £(Si) - E(TT) [eV] 



FIG. 6. The same as in Fig. [3] but for population dynamics calculated by 
the NIBA-type master equation, which is perturbative to second order in the 
electronic couplings, Vij. 



see that the rate of fission decreases drastically as r] — » 0, ap- 
proaching a 100 ps timescale for 77 = 0.2. The fission rates 
(obtained by a numerical fit of the TT population growth) are 
then plotted in log-log scale, Fig. 13b), very clearly confirm- 
ing the superexchange scaling k oc rf . Another interesting 
feature is apparent in the long time dynamics of Fig. Ha). 
To quantify this behavior, in Fig. EJc), we plot the equilib- 
rium population of TT, obtained as PttO — » 00) (equivalently 
Z~i(TT\exp(-H e ilkBT)\TT)) as a function of the electronic 
coupling strength rj. Clearly, for increasing 77, the equilibrium 
population shows a noticeable decline, which is straightfor- 
wardly explained: for stronger values of coupling, the zeroth- 
order diabatic states are more strongly mixed, such that the 
lowest energy adiabatic eigenstate develops a larger fraction 
of non-7T states, effectively depleting the diabatic TT popu- 
lation. 

The mechanistic features evinced in this and the previous 
section constitute the main results of this paper. Namely, a su- 
perexchange two-electron-transfer phenomenon, utilizing vir- 
tual states in the CT vibronic manifold, is entirely consis- 
tent with the observed features of singlet fission in pentacene. 
As such, a CT-mediated singlet fission mechanism cannot be 
ruled out based solely on the argument that CT states are too 
high in energy. From a more technical point of view, only 
a dynamical master equation which can account for fourth- 
order effects in the electronic couplings, such as Redfield the- 
ory, is able to correctly predict this behavior. Before consid- 
ering the feasibility of a direct mechanism, mediated by a true 
two-electron coupling matrix element, we investigate in the 
next section the role played by the phonon bath degrees of 
freedom. 




0.1 1 0.2 0.4 0.6 0.8 1.0 

Electronic coupling strength, r\ Electronic coupling strength, 77 

FIG. 7. Dramatic slowing down of singlet fission dynamics for decreasing 
electronic coupling strength 77 (a); note that the time axis is in log-scale. The 
numerically extracted fission rate obeys the predicted superexchange scaling 
k ~ if (b), however the equilibrium population of TT decreases with increas- 
ing coupling, due to enhanced mixing with non- TT states (c). 



C. The effect of the bath 

Within the weak system-bath coupling approximation of 
Redfield theory, the adiabatic population transfer rates k a ^>p = 
Rpp aa are given by 

k^ = h~ l C a/3 J(o) a/3 )n(a) a/3 ) (14a) 
k^ a = h~ l C a pJ(cL> a p) \n(u> a p) + l] (14b) 

where C a p = £i |(/3|/)| 2 |(tf|/)| 2 arises from the change of basis, 
n(co) - [exp(tiCL)/kBT) - 1]~ is the Bose-Einstein distribution, 
and we assume cj a p = cj a -cjp > 0. Clearly these transfer rates 
satisfy the detailed balance condition, 

-—^ = exp (tiCL> a plkBT\ . (15) 

kj3^>a 

Furthermore, the transfer rates are directly proportional to the 
spectral density evaluated at the eigenvalue energy difference. 
Physically, the transfer rate between two states depends on the 
availability of strongly-coupled phonon modes at the required 
energy difference as well as their thermal occupancy, so that 
phonon absorption and emission facilitates the electronic en- 
ergy transfer. 

The fission rate will clearly depend on the parameters and 
functional form of the spectral density, J(co). As discussed 
previously, the chosen reorganization energy A = 50 meV is 
conservatively small, and so here we analyze the dependence 
of the fission rate on the strength of the system-bath coupling, 
quantified by A. Because J(oj) oc X for any form of the spectral 
density, Redfield theory predicts the trivial linear dependence 
^fiss 00 ^ via Eqs. (IT4b : see Fig. [8j While this is indisputably 
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50 100 150 200 
Reorganization energy X [meV] 

FIG. 8. Calculated fission rate for a pentacene dimer with varying system- 
bath coupling, quantified by the reorganization energy, A. Secular, Markovian 
Redfield theory (filled circles) predicts a linear dependence, which is known 
to be accurate for small A but becoming more inaccurate for large A (indicated 
by the shaded region). Realistic values for pentacene are A « 50-150 meV, 
which reliably predicts a fission rate k « 2-10 ps _1 i.e. r « 100-500 fs, in 
reasonably good agreement with experimental rates of 80-200 fs. 



the correct behavior in the small A limit, it becomes incor- 
rect for A sufficiently large. 35 This breakdown is visually sug- 
gested in Fig.[8l with a shaded region indicating at what point 
the theory may become inaccurate, initially quantitatively but 
ultimately qualitatively 46 

We conclude this section by considering different func- 
tional forms of the spectral density, although we stress that 
the Redfield formalism is equally accurate and efficient for 
any form of the spectral density, including ones determined 
numerically. In addition to the Ohmic form utilized so far, we 
will also consider a super-Ohmic spectral density of Debye- 





50 100 150 200 250 300 350 400 

300 
- 240 > 
180 B 
120 § 

60 il| 


50 100 150 200 250 300 350 400 
Tiuj [meV] 

FIG. 9. Three different forms of the spectral density investigated here, 
along with the electronic eigenvalue differences for pentacene (orange ver- 
tical sticks). The overlap between these energy differences and the spectral 
density, i.e. the ability to absorb and emit resonant phonons, largely deter- 
mines the rate of population transfer and hence singlet fission. 



type phonons with exponential cutoff, 

Jffico) = ^^^ 3 exp(-3o;/D i7 ), 

ij 



(16) 



with A = 50 meV and Q = 180 meV as before. Lastly, we will 
consider a broadened stick- spectrum of pseudo-local phonon 
modest 



1 E- 



AkjjFkjjOJ 



2 ' 
k,ij 



(17) 



with total reorganization energy A t j = Yuk ^k,ij- The latter has 
three phonon modes, chosen based on the results presented 
by Girlando et al.^ with A k = {15, 20, 15} meV (£ k A k = 50 
meV) at frequencies ficok = {150, 175, 200} meV, with uniform 
broadening HTk = 50 meV. 

As just discussed, the rate is largely determined by the over- 
lap of the spectral density with the eigenvalue energy differ- 
ences. In Fig. [3 we plot these three different spectral den- 
sities along with the energy differences, visually portraying 
which modes mediate which electronic transitions. Clearly, 
the overlap is most uniform for the Ohmic and least uniform 
for the pseudo-local spectral density: while the pseudo-local 
spectrum greatly enhances some electronic transition rates, it 
greatly diminishes other, in particular those with relatively 
small energy differences. In Fig. [9] we also show a possible 
Peierls-type (off-diagonal) spectral density based on the crys- 
tal phase calculations of Girlando et al.r^ though we reiterate 
that such coupling was not employed in the results presented 
here. Based on the overlap argument given above, it is vi- 
sually apparent why these low frequency fluctuations do not 
efficiently mediate the fission process, as we have found in 
dynamics calculations not shown. 

The fission dynamics due to these various spectral densi- 
ties are found to be only mildly different, as shown in Fig.fTOl 
This behavior is due largely to the uniformly strong coupling 
to phonon modes near 160 meV as these phonon modes me- 
diate the important S i => (S i ^ CT) and (S i ^ CT) => TT 
transitions, see Fig.H Based on these results, one can in prin- 
ciple imagine devising clever ways to engineer either the elec- 
tronic spectrum or the phonon spectrum to realize maximal 
singlet fission rates. In practice, an accurate determination of 
the "correct" spectral density is very difficult, but because it 
can clearly shape the observed electronic dynamics, we con- 
sider this an important topic worthy of further study. 



D. Direct mechanism 

While the results so far have focused on the features of a 
CT-mediated singlet fission mechanism, we now briefly inves- 
tigate the feasibility of a direct-coupling pathway. Ultimately, 
we find that the required two-electron integral which couples 
S\ to TT is too small to explain efficient singlet fission in 
pentacene dimers Of course, it remains to repeat the analysis 
in crystals, for which the multiexcitonic density of states is 
enhanced as compared to the dimer. This calculation will be 
reported in a future paper. 
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FIG. 10. Singlet fission population dynamics for the three spectral densities 
depicted in Fig. [9] 



Using the methodology described in Sec. HH we calculated 
the direct coupling matrix element, 

(S 1 S \6ei\TT)= ^{(L A L B \H B L A )-(H A H B \L B H A )}, (18) 
where 

(ij\kl) = J J 3 n J d 3 r 2 0;(ri)05(r 2 )ri 2 Vit(ri)0/(r2), (19) 

finding it to be less than 1 meV. This value (which is in agree- 
ment with more accurate calculations using block-localized 
DFT 49,50 ) should be contrasted with the effective superex- 
change coupling, Eq. dTQb . which is approximately 50 meV 
for the pentacene dimer. 

These qualitative arguments are numerically confirmed in 
Fig-El where we show the singlet fission dynamics in the ab- 
sence of a CT-mediated pathway, i.e. the electronic coupling 
to all CT states is zero, thereby reducing the number of rele- 
vant states to three. In particular, we use E(S \)-E(TT) = 250 
meV, A = 50 meV, and bath frequency and temperature as be- 
fore. By varying the strength of the two-electron coupling, 
we observe that the direct mechanism only becomes compet- 
itive at unrealistically large values. As such, we conclude 
that, while feasible, the direct singlet fission mechanism is 
subdominant as compared to CT-mediated superexchange, at 
least for pentacene dimers. 



E. Covalently linked dimer 

In this final section, we briefly investigate a molecular ge- 
ometry that is different from those considered before. In par- 
ticular, inspired by the synthesized tetracene molecules of 
Miiller et al.riLSI we consider the analogous pentacene dimer 
shown in Fig. [121 In Refs. [13 and 0, the authors found 
that less than 1% of singlets underwent fission, which is in 
stark contrast to the relatively high fission yield of crystalline 
tetracene. 

The covalently bound dimer pair contains a benzene linker, 
which can facilitate a through-bond coupling in addition to 
the through-space coupling available in the crystal. We 
can straightforwardly estimate the through- space coupling 



1.0 




1 2 3 4 5 

Time t [ps] 

FIG. 11. Singlet fission population dynamics in the absence of CT states, for 
varying values of the direct electronic coupling element given in the legend. 
Sub-picosecond fission is only observed for the unphysically large value of 10 
meV, to be contrasted with theoretical estimates ranging from 5 to less than 1 
meV. 



with the methodology described above by excising the ben- 
zene linker and terminating the dangling bonds with hydro- 
gens. This calculation yields one-electron coupling matrix 
elements, ty, of approximately 0.02 meV, which is three or- 
ders of magnitude smaller than in the native crystal confor- 
mation. The mixed through-bond couplings, t HL and t LH , are 
not so easily obtained, but the HOMO-HOMO and LUMO- 
LUMO couplings can be estimated using the energy- splitting 
method alluded to previously. For this calculation, the ben- 
zene linker is retained in the molecule, and the splitting of 
the HOMO and LUMO orbitals is assigned to twice the re- 
spective coupling. We calculate this through-bond coupling 
to be on the order of 20 meV, significantly larger than the 
through-space contribution, but still 4-10 times smaller than 
the through-space coupling of the crystalline dimer pairs in- 
vestigated above. In light of the scaling analysis presented in 
Sec. IIIIB1 it is clear that such a reduction will yield a signif- 
icantly slower rate of fission, increasing the timescale from 
sub-picosecond to as much as one nanosecond, extending the 
timescale over which fission must compete with other decay 
mechanisms. The torsional modes about the benzene linker, 
which couple to these off-diagonal electronic hopping matrix 
elements, are expected to be of relatively low frequency and 
thus would not qualitatively affect the fission dynamics, simi- 




FIG. 12. Two different views of a covalently linked pentacene dimer akin to 
the tetracene dimers of Refs. l20l and [2]]. The electronic couplings are esti- 
mated to be significantly smaller than in the crystal phase dimer pairs inves- 
tigated above, predicting a much smaller rate of fission. 
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lar to the analogous observation we made for the crystal. This 
example calculation provides a possible explanation for the 
low observed fission yield 20,21 and encourages efforts to ex- 
plore electronic coupling effects through a combination of co- 
valent bonding and optimization of geometric orientation. 



IV. CONCLUSIONS 

In this and the previous work, we have presented and ap- 
plied a unified microscopic theoretical framework for the in- 
vestigation of singlet fission electronic structure and dynam- 
ics. We have emphasized the role played by molecular vi- 
brations or phonons in mediating non-resonant excited state 
energy transfer. Such finite temperature relaxation mecha- 
nisms are responsible for population transfer, coherence de- 
phasing, and eventual thermalization. These processes may be 
numerically studied and quantified using a perturbative quan- 
tum master equation subject to its regime of validity. We have 
justified and pursued Redfield theory for singlet fission chro- 
mophore systems, advocating its favorable trade-off between 
accuracy and efficiency. Within our framework, we have per- 
formed a thorough investigation of singlet fission in pentacene 
dimer systems, in particular investigating the real-time fission 
dynamics. In principle, all parameters needed to carry out the 
real-time quantum dynamics at the level detailed in this work 
may be estimated from microscopic considerations. Instead, 
in this work we have carried out a systematic variation of sev- 
eral parameters which influence singlet fission efficiency, pro- 
viding insight into the underlying mechanistic details which 
are beyond the reach of experiment. 

Most notably, we have provided evidence for a CT- 
mediated superexchange mechanism in pentacene dimers, 
which is more efficient than a direct mechanism, even in 
the presence of very high-energy intermediate CT states. 
These results are in stark contrast with previous theoreti- 
cal predictions. 23,24 However, our superexchange mechanism 
bears many similarities to the conventional direct mechanism 
and it is in some sense a hybrid between the two competing 
mechanisms. We have investigated the way fission rates and 
yields are modified by shifting the electronic energy levels and 
scaling the electronic couplings. Our study of the role of the 
phonon bath properties underscores the importance of having 
resonant phonon frequencies to mediate efficient energy trans- 
fer. Lastly, we have shown that, at least in pentacene dimers, 
the direct coupling singlet fission pathway is extremely inef- 
ficient and would require an unphysically large two-electron 
matrix element to compete with the mediated superexchange 
mechanism whose timescale is already in very good agree- 
ment with experimentally observed fission rates in pentacene 
crystals. 

It should be emphasized, however, that it may be dangerous 
to draw conclusions about singlet fission in the bulk based on 
the calculations for idealized dimers presented here. While 
our work clearly shows that the existence of CT states that 
are energetically high-lying in no way obviates their impor- 
tance in singlet fission, it does not prove that superexchange 
is relevant for bulk pentacene or other singlet fission mate- 



rials. Further, the possibility exists that an enhancement of 
the direct mechanism may take place in bulk materials due 
to the augmented density of multi-exciton states expected in 
large clusters and bulk solids and films. Our preliminary cal- 
culations on larger pentacene clusters based on the formalism 
presented here, which allows for the efficient simulation of 
hundreds of quantum states, suggest that many of our conclu- 
sions are unaltered, in particular that CT states are intrinsically 
connected to fission and the predominant mechanism is not a 
direct two-electron process. This work will be presented in a 
forthcoming paper. 
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